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Abstract 

We present the results for the fuh electroweak two-loop bosonic contributions to the 
effective leptonic mixing angle of the Z boson, sin^^eff; in the Standard Model. A 
method applied to extract collinear divergences from two-loop vertex-functions is 
described. Comparisons of our results with those from a recent previous calculation 
show complete agreement. 



1 Introduction 

Measurements of the various asymmetries of the Z resonance have determined the effec- 
tive leptonic mixing angle sin^^efr with high accuracy. The current experimental value 
is 0.23153 ± 0.00016 [1]; a linear electron-positron collider with GigaZ capabilities could 
even reach an accuracy of 1.3 x 10~^ [21 [3]. Comparison with the theoretical prediction 
of the Standard Model yields stringent bounds on the Higgs-boson mass Mu- The impor- 
tance and precision of sin^6'efr requires sufficient control on the theoretical accuracy of the 
prediction by providing adequate higher-order calculations. 

sin^6'eff is determined from the ratio of the dressed vector and axial vector couplings 
gy^A of the Z boson to leptons [4j, 

sin^^efT = ^(l-Re^). (1) 

In the on-shell renormalization scheme, it is related to the vector-boson mass ratio or, 
equivalently, to the on-shell quantity = 1 — M^/M| via 

sin^^off = ^^s^ ■, k = 1 + Ak , (2) 



involving the k factor, which is unity at the tree level and accommodates the higher-order 
contributions in Ak. My^ can be related to the precisely known Fermi constant via the 
relation 

2 A vra 



where Ar summarizes the higher-order contributions. Beyond the one-loop order, the 
universal QCD corrections [5|[6], the complete electroweak fermionic [7] and bosonic [8] two- 
loop corrections, as well as leading higher-order contributions via the p-parameter [9l [TO] 
and the S'-parameter [11] are known for the Standard Model. 

The universal higher-order terms to the self-energies from QCD and through the p- 
and S'-parameter enter also the quantity Ak in ([2]). The fermionic electroweak two-loop 
corrections (involving at least one closed fermion loop) have been established [T2l [T3] ; the 
M//- dependence of the bosonic corrections to Ak was given in [H], and only recently a 
complete calculation of the bosonic contributions was reported [15] . 

In this paper we complete our result for the bosonic two-loop contributions by ac- 
complishing the M/^-independent part of Ak, which was missing in [T3]. Together with 
our result from [Il|, this yields an independent result for the complete bosonic two-loop 
corrections for sin^6'eff, based on different methods. 
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2 Calculational strategy 



Expanding the dressed couphngs in ([T]) according to gv,A = Qva \ ^ + 9va + 9va + ' ' " j iii 
powers of a, the O (a^) contribution to sin^^efr in the loop expansion 

sin^^eff = sm'e^^ + + sm'e^^ + O {a') (4) 

is obtained, 




Hence, besides two-loop diagrams, also products of one-loop contributions have to be taken 
into account. They play an important role in the cancellation of IR-divergences (see sec- 
tion [3]). But the major task consists in the calculation of irreducible two-loop Zii-vertex 
diagrams {i = e, to be precise). The Z-boson couplings in ([T]) appear in the renormalized 
Zi£ vertex for on-shell Z bosons; for the two-loop contributions entering the expression ([5]) 
we need the renormalized two-loop vertex 

ff^(^)(M|) = 7. {g'^'-g'?!.) ■ (6) 

As done in [T3] we split the renormalized vertex into two UV-fimte pieces, 

= K^^^^o) + rf ] + K"(^)(M|) - rf ^(^)(o)] . (7) 

pZ«(2) denotes the corresponding unrenormalized Z£i vertex for on-shell leptons and 

momentum transfer = M|, and T^'^ is the two-loop counter term. Details on the 
renormalization are given in [13], where the expressions are general and comprise also the 
bosonic two-loop contributions. 

The first term in ([7j) can be computed as in the fermionic case, which means generating 
Feynman diagrams with the help of FeynArts [TB] and applying TwoCalc [T7] to reduce 
the amplitudes to standard integrals. The resulting vacuum integrals are calculated using 
analytic results |T8til9j, whereas the two- loop self-energies with non- vanishing external mo- 
menta, which appear in the counter terms, are calculated with the help of one-dimensional 
integral representations [20j or the methods described in [21] . 

The Ji?-finite contributions to the second term in ([7]) can be calculated in analogy to 
the Mff-dependent part of the bosonic corrections [H]. This means that vertex corrections 
are computed applying the methods from ^21j together with further improvements in order 
to increase the numerical stability. Non-planar diagrams are calculated using the method 
described in [15], and diagrams with self-energy insertions are computed using dispersion 
relations as described in [T3]. 
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(a) (b) 




(c) (d) 
Figure 1: Diagrams with IR-divergences 

3 Treatment of IR- and collinear divergences 

As a new feature compared to the M/^-dependent subset of the bosonic corrections, dia- 
grams with internal photons appear that can be IR-divergent. These divergences, however, 
cancel in the complete result rendering sin^^eff as a IR- finite quantity. 

The IR-divergent diagrams are shown in Fig. [H The grey circles represent one-loop 
insertions. In order to verify the cancellation of these divergences, we have used the 
methods described in [22]. These methods allow to extract all appearing IR-divergences 
in terms of IR-divergent one-loop integrals. After this extraction the cancellation of these 
integrals can be checked analytically. The cancellation occurs between the various two-loop 
diagrams within the set of Fig. [I]^a)-(c), and between the two-loop diagrams of Fig. [U^d) 
and the product of one-loop diagrams occuring in This is shown schematically in 
Fig.H 

Moreover, collinear divergences appear. In order to regularize them we have kept 
the electron mass me finite where necessary. In addition we have expanded the resulting 
expressions in me, so that the divergent behavior shows up as terms proportional to In^(mg) 
and In(mg). Afterwards it was checked analytically that the terms proportional to In^(mg) 
cancel, whereas the cancellation of the ln(mg)-terms was checked numerically. Our methods 
to extract these logarithms are explained in the following section H] by means of the two 
non-planar diagrams shown in Fig. [3l 
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+ finite terms. 



Figure 2: Compensation of IR-divergences 



4 Extraction of collinear divergences 

CoUinear divergences can show up in massive diagrams which do not contain any soft 
divergences, or they can overlap with soft divergences. In both cases we have managed to 
extract exphcitly the coefficients of the logarithms in (which is the natural regulator for 
these divergences) and the constant term. For overlapping divergences we have used the 
results of [22j, while for the other ones we have used a subtraction method in parametric 
space. In order to explain it we consider the two collinear configurations of the non-planar 
diagrams in Fig. [3], which arise in the computation. 



Figure 3: Examples for coUinearly divergent diagrams. All momenta are incoming, and 
P = P1+P2. 

In momentum space the generic diagram contains the following scalar and tensor integrals. 



-P 




222 
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[1] = ql - ml, [2] = (gi - p^Y - ml [3] = {q^ - q, + p.f - ml 
[4] = (gi - g2 - P2)^ - [h\=ql-ml [Q] = [q^ - p^f - ml 



Following the discussion given in [21], we first combine the propagators [1] and [2] with 
a Feynman-parameter zi, the propagators [3] and [4] with a Feynman-parameter Z2-, and 
the propagators [5] and [6] with a Feynman-parameter z^. Then we combine the qi and 
5i — Q2 propagators with a parameter x. After the gi-integration we combine the residual 
propagators with a parameter y and carry out the 52 integration. For the scalar integral 
V^22 5 the resulting expression reads (we have changed x ^ 1 — x with respect to [2T]): 

^222 = j dxdy dzi dz2 dz^ xy{l — x){l — y) f/^^, (9) 
Jo 

U = xy{l-x){l-y)i + y{l-x)x{,zi]pl]mi,m2) 

+ xyx{z2\P'^\m^,m^) + x (l-x)(l -y) ^(zs; p^; mg, mg), 

C = - [i^^ (22-2:3) {I-Z1-Z2) +pI {z2-Z^) (I-Z1-Z3) {l-Zi-Z^) {l-Zi-Z2)\ , 

x{z]p^]m,M) = -p^ z{l-z) {^-z) + M"^ z . 



Here and in the following we use the short-hand notation 

1 pi i>i i>i 

dxdy ■ ■ ■ dzs = dx dy ■ ■ ■ dzs . 

Jo Jo Jo 



4.1 Non-planar diagram: two internal fermions 

In Fig. [3](a) the following mass configuration of the V222-faniily appears, 

mi = ms = 0, 1712 = ^e, ^3 = = ^6 = = M|, vi = vi = ''^1 ■ (10) 

The divergence comes entirely from the fermion-photon interaction, so it is possible to set 
p\ = Q. Inserting these values into ([9]) the resulting expression for the scalar integral ¥2^ 
reads 

(1) /"^ 2^(1 - ^) 
^22 = I dy dzi dz2 dzs y {I - y) /, with I = dx — ^ , (11) 

Jo Jo [x [ax + 0) c{x)\ 

where a, b and c(x) depend on Mw, P^, and on the Feynman-parameters y, zi, Z2, and 
Zs- In particular, b is linear in zi and z^, 

b = aziz^ + Pzi + 72:3 + 5, with 

a = Mly{l-y), 

f3 = -Mly{l~y)z2, 

7 = {l-y)[M^-Mly{l-Z2)], 

S = y[M^-Mlyz2{l-Z2)]. (12) 
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Moreover, the following relations hold, 

c{{))=yzl a + b = y[M^-Mlz2{l-Z2)]. 



(13) 



From ( fTTl) one can see that, for mg = 0, J is divergent at x = 0. The divergence can be 
extracted as In ml in the following way: 



la + h 



1 

dx 







x(l — x) 



X 



[x {ax + b) + ml c(a;)]^ [xh + ml c 



[xb + ml c{0)Y 



(14) 



la is not divergent and we are free to set m^ = 0, 



dx — 



X 



1 — X 

[ax + by 



1 

62 



1 /, a + b 



1 . 



In lb, the divergence can be extracted by integrating in x, 

1 



h = - 



62 



Inmg + Inc(O) - ln6 + 1 



(15) 



(16) 



Thanks to the linearity of b with respect to zi and z^, we can now perform the zi,Z3 
integrations explicitly. Taking into account ( fT2l) and (|T3l) . the following types of integrals 
remain: 



I 



dzi dz3 — 
" 



(i2;i (i2;3 



iln|l + 



d 



d 



Li2 



[75] 



Li2 







' , ln6-l 
d^i dz3 — - — 
o2 



- ln[a/?7(5] In (l 



Lig 



— Lig 



d 



a [yd] 



^^)+Hm Infl 



+ In [75] In 1 



d 



a [y6] 



-ln6 In (1-4 

ao 



Li, I — 

a 



(17) 



where we have introduced d = a6 — P'j and [ai...Q;„] = ai + ... + an- It is important to 
note that the only denominator appearing in the result {d) is now multiphed by "regulator 
functions" which vanish when d goes to 0. The result is therefore smooth enough to be 
directly inserted into fllip . and the remaining Z2,y integrations are done numerically. 

The computation of the tensor integrals V^g^''^'^'^'^''^ is processed in the same way. The 
new feature is the presence of z^ to some power m > in the integrand. The previous 
formulae are now replaced by the following ones, where an integration by parts in z^ is 
performed: 
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b' Jo d V bS] {[a-f]zs + mi 



\z dz z-^""'' - r dz "^'^'^ lli ( ["^n Li ( ""^ + ^ 
dz,dz,z,^ - ^^3^^|Li,^-j-^J-Li,^-— ^ 







IP' Jo d \ \a[al3'-)S\J \a{[a'-)]z^ + [j38\) 

Li2 { -4^ + Li2 { ^-^-ttI - ln[a/575] In fl • 



+ ln([a7]2;3 + [/55]) In (l -. — / .o^^, 

+ ln[7(5] In fl - ln(723 +5) In fl ■ 



a;[7(5]y \ a (7^3 + 5)^ 

4.2 Non-planar diagram: four internal fermions 

Fig. [3](b) corresponds to the configuration 

c\ c\ ry ry O / \ 

mi = 0, 1712 = rris = = mg = me, = Mz P = M^, Pi= P2 — ^e- (19) 

This is an example for a diagram with a double collinear divergence. In order to regularize 
these divergences the masses of all internal electrons are kept and p\ = = is used. In 
the following the scalar integral V222 is considered; the tensor integrals can be calculated 
as in section 14.11 . 

We apply the same parametrization as in section 14.11 and change parameters according 
to 2:2,3 — 1 — Z2,3, z\ zz- The following expression is obtained, 

V222 = / dydzsy{l-y) J, with J= dxdzidz2- — ■ ^ r^, (20) 

Jo Jo [p + m^a — 10) 

P = -X {yz2{l - Z2) - {l-x){l-y)[zi + y{zi-Z2){z2-Z3)]} Ml , 

a = -xy {l-x){l-y){zi-Z3f + {l-x)yzl + xy + X {l-x){l-y) {1-zif . 

Since p is of the type p = x {zi A + Z2 B + zi Z2 C) , one can easily see from eq.(l20l) that, 
for me = 0, V222 diverges at x = and at zi = = Z2. J can be calculated in the following 
way: 

J = Ja + Jb + Jc + Jd = [J — Jx — Jz + Jxz] + [Jx — Jxz] + [J z — Jxz] + Jxz , (21) 
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dx dzi dz2 



xil — x) 



{p + ml ay 



J. 



dx dzi dz2 



X (1 — x) 



zf=Z2=ziZ2=0 



J X 



dx dzi dz2 



X (1 — x) 



(p + ml a) 



x'^=zf=Z2=ziZ2=0 
xm'^=zi m^=Z'2 m^=0 



J a is finite for m^. = 0, while Jf, and Jc involve a single divergence (at x = and zi = = Z2, 
respectively), giving raise to a simple logarithm in m^. Finally, Jd contains the double 
divergence at x = and zi = = Z2, and we extract also a In^ m^ term accordingly. Now 
we define: 



a + b = -y Z2{l-Z2)Ml - i5, 

b = {{l-y)[zi + y{zi-Z2){z2-Z3)]-yz2{l~Z2)]}M^ 

ao + bo = -y Z2 M| - i6, 

bo = {i^-y)[zi -yz3 izi-Z2)] -yz2]Ml - i6, 



(22) 



i6. 



A=il-x)il-y)il-yzs)Ml-t5, 
Ao = A\,=o = {l-y){l-yz3)Ml - 



i6. 



B = -y[l- (l-x){l-y)zs\M^ -i6, 



=0= [l-x{l-y)] [x+{l-x)yzl], 

- O'\x=zi=Z2=0 

With these notations we obtain for J^, setting m^ 



a 



a 



Zl=Z2 = 



x=Q 



y4 



0, 



Ja 



dx dzi dz2 



X 



X 



{ax + by 



1 



1 — X 



+ 



1 



dzi dz2 



In 6 - ln(a + b) - I 
62 



{aQX + b^y ^0 

ln6o - ln(ao + 6o 
1? 



(23) 



The resulting expressions are of the form flTTl) and can be integrated in zi and z^ accordingly. 
Since in this case both 7, 5 and d of (ITTl) factorize a Z2, we get a factor l/z2, that is not 
multiplied by regulator functions. However, even if the integrals containing 1/6^ and l/fcg 
are separately divergent, the combination of the two is finite, because they have the same 
behaviour around Z2 = 0. Therefore, after applying (fT7|) . the contribution from J a to V2II 
can be safely integrated numerically. 

For Jfo we proceed in an analogous way integrating in x, and obtain 



dx dzi dz- 



X 



X 



dzi dz2 



Xxb + mlaoof (xbo + mlaoof _ 
In 6 — In (Too ^1 In 60 — In cxoo — 1 



62 



— ln(m^ 



1 

62 



62 



(24) 
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Again, these terms are of the same form as in the expression for Ja and can be integrated 
analytically in zi and 2:3 and then numerically in y and 2:2- 

To extract the logarithmic behaviour for Jc we have to integrate in zi and 2:2, yielding 



J. 



dx dzi dz2 



dx — 

X 



x(l — x) 



X 



[x{Azi + Bz2) + mlaoY [x{AqZi + ^02:2) + mlaoof 
l-x) ^n{A+B) - In A - In 5 - In X + In ctq] 



AB 



hi{AQ + Bq) — In y4o — In _Bo — Ina; + Incroo 
A^^ 



+ ln(mg 



1 — X 



1 



AB AoBo 



(25) 



The factors y, {1 — y), and (1 — x) in the denominators of this expression are always 
cancelled by those present in the numerator. The overall factor 1/a; is again regularized by 
the difference between the terms with A, B, o"o and those with Aq, Bq, aoo- The remaining 
y and z^ integrations can be done numerically. 

Finally, for we integrate in x, zi and 2:2, yielding 



J, 



dx dzi dz2 



Li2 



[x{AoZi + ^02:2) + mlaooY 



AoB, 







Ao + B^ 



Lio 



An 



The dilogarithms can be expanded according to 

^2 1 



Li2 - 



c 



TC 

'~6 



Li2 



5, 







m^e^^OQ 



2 {\nc-\nmiy + 0{mi), 



(26) 



(27) 



giving rise to a term proportional to In m^. The resulting expression is suitable for nu- 
merical integration. 



5 Results 

In Tab. [1] the set of input parameters for our numerical evaluation is listed. Mw and 
Mz are the experimental values of the W- and Z-boson masses, which are the on-shell 
masses. They have to be converted to the values in the pole mass scheme [7], labeled as 
Mw and Mz, which are used internally for the calculation. These quantities are related 
via Mw,z = Mw,z + ^wz/C^ Mw,z)- For Tz the experimental value (Tab. [I]) and for Tw 
the theoretical value has been used, i.e. T]y = 3 G^M^r/ (2-\/27r) (1 + 2as (M^) / (Svr)) 
with sufficient accuracy. 

Our results in terms of An for the bosonic contributions are shown in Tab. [2], in com- 
parison with the fermionic ones. The last two columns of Tab. [2] contain the corresponding 
results from [T5], which agree with ours. 
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parameter value 
3% 80.404 GeV 

Mz 91.1876 GeV 

Tz 2.4952 GeV 

mt 172.5 GeV 

Aa (M|) 0.05907 
a, (M|) 0.119 

1.16637 X 10-5 

Mv^ 80.3766 GeV 

Mz 91.1535 GeV 

Table 1: Input parameters. Mw and Mz are the experimental values for the W- and Z-boson 
mass, whereas Mw and Mz are calculated quantities in the pole mass scheme. 



Mh [GeV] 


A4r2 X 10- 


A-s:^ X 10- 


AntZ X 10- HZ 


a4::) X 10- [15] 


100 


1.07 


-0.74 


1.07 


-0.74 


200 


-0.33 


-0.47 


-0.32 


-0.47 


600 


-2.89 


0.18 


-2.89 


0.17 


1000 


-2.62 


1.11 


-2.61 


1.11 



Table 2: Bosonic two-loop corrections to Ak in comparison with the fermionic ones, and the 
results from |15j . 



One can see from Tab. [2] that the bosonic corrections are basically of the same order of 
magnitude as the fermionic ones, as long as the IV-mass is taken as a fixed input parameter. 
The situation changes, however, when instead of the l^-mass the Fermi constant is taken 
as an input parameter and M]^ is derived via ([3]). In order to identify the various sources 
of the two-loop contributions to sin^^cfr, "we expand both quantities Mw and k according 
to {Mw, /t} = {Mw, + {Mw, + C(a^), yielding 

sin^^eff = sin2^fg'''+"^ + A sin^^efflAM^^'^) + A sin^^efrlAfi:^"')) + 0{a^) (28) 

1 - , > «:(*-^+-) - 2 ^ AmI:'^ + ( 1 - ^ ) A«:('^^) + 0{a') 

M| J ml \ ml) 

The first term represents the one- loop result, while the other two terms correspond to 
the electroweak two-loop contributions to the shift in the W mass, AM^ \ and to the 
shift in K, Ak^" \ As can be seen in Tab. |3l the fermionic corrections from M^/ dominate, 
while the contributions from My/ and k in the bosonic sector cancel to a large extent. 

In Tab. |U the results for M^y, as taken from [TJ], and sin^^eff are shown with only the 
fermionic two-loop corrections (M^^""™ and sin^6'gg™) and with both fermionic and bosonic 
two-loop corrections (M^™+^°' and sin2^||g™+^°') for various values of Mr. Besides the 
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Mh 


Asin2^eff(AM(^;L) 




Asin2^eff(AM(^l) 


Asin^^^efrlA^::)) 


100 


93.89 


2.38 


1.93 


-1.65 


200 


98.51 


-0.73 


0.97 


-1.05 


600 


106.89 


-6.43 


0.19 


0.40 


1000 


105.36 


-5.83 


-1.16 


2.47 



Table 3: Fermionic and bosonic electroweak two-loop corrections to sin^^cfr- Mh is given 
in GeV; the values for Asin^6'efr have to be multiplied by 10~^. 



Mh [GeV] 


M^™ [GeV] 


^term+bos jQ^yj 


sin^^ff"^ 


sin^Cf?"^'" 


sin^^^f [10-^] 


100 


80.3694 


80.3684 


0.231459 


0.231461 


0.02 


200 


80.3276 


80.3270 


0.231792 


0.231792 





600 


80.2491 


80.2490 


0.232346 


0.232352 


0.06 


1000 


80.2134 


80.2141 


0.232587 


0.232599 


0.12 



Table 4: My/ and sin^^cfT without (ferm) and with (ferm+bos) bosonic two-loop corretions. 
s\v?9^^^ is the purely bosonic contribution to sin^^efj . 



electroweak one- and two-loop corrections, the results for sin^^eff also contain the QCD- 
corrections of order 0{aas) and O^aaj.) p], as well as the leading three-loop corrections 
of order 0{a'^asmj) and 0{a^mf) [9]. 

The shifts in sin^^efr originating from the full set of bosonic two-loop contributions are 
displayed in the last column of Tab. HI According to the cancellations mentioned above, 
they are much smaller than the corresponding shifts for a fixed value of the W mass as an 
input parameter, which would result from eq. ([2]) with the values of An given in Tab. [2l 

6 Conclusion 

In conclusion, we have completed the calculation of the electroweak bosonic two-loop cor- 
rections to the effective leptonic mixing angle sin^6'efr. A comparison of our results with 
the ones in [15] was performed and full agreement was found. The bosonic corrections are 
basically of the same order of magnitude as the fermionic ones when the ly-mass is taken 
as a fixed input parameter. If instead the Fermi constant is taken as input, with as a 
derived quantity, the bosonic corrections to sin^^gff are quite small. 

In addition, the use of the general techniques of [22] to treat IR-divergences provides a 
good test of their applicability to physical processes. In particular we have shown that in 
this approach the cancellation of the infrared poles can be verified analytically. Moreover, 
our complete calculation of the two-loop electroweak corrections of sin^^es is the first 
physical application of the numerical methods of [21] showing that these methods are 
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really applicable in practice for the calculation of physical quantities. In particular within 
this approach massive two-loop vertices with many mass scales are computed "directly" , so 
not using mass expansions, which could represent a good way to deal also with calculations 
beyond the Standard Model, where often many (unknown) massive particles come into the 
game. 
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